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ABSTRACT 



We used HARPS to measure oscillations in the low-mass star r Cet. Although the data were compromised by instrumental noise, we 
have been able to extract the main features of the oscillations. We found r Cet to oscillate with an amplitude that is about half that of 
the Sun, and with a mode lifetime that is slightly shorter than solar The large frequency separation is 169 /jHz, and we have identified 
modes with degrees 0, 1, 2, and 3. We used the frequencies to estimate the mean density of the star to an accuracy of 0.45% which, 
combined with the interferometric radius, gives a mass of 0.783 ± 0.012 Mq (1.6%). 
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1. Introduction 

In the past few years, a new generation of high-resolution, high- 
precision spectrographs has been providing unprecedented op- 
portunities for studying the fine details of stellar interiors and 
evolution through the detection of tiny stellar oscillations. The 
observation and analysis of stellar oscillations, or asteroseismol- 
ogy, has the potential to change dramatically our views of stars. 

The G8 V star t Ceti (HR509; HD 10700; HIP 8102; V = 
3.50) is expected to have a convective envelope and therefore to 
display solar-like oscillations. Since t Cet has a lower metallic - 
ity than the Sun ([Fe/H] = -0.5 + 0.03: ISoubiran et al.llI998h . 
it bridges the gap towards very metal-poor population II astero- 
seismic target stars s uch as v Ind, where solar-type oscillation s 
have been detected jBedding et all 120061; ICarrier et all |2007|) . 
Moreover, among stars for which a dete ction of solar-type os- 
cill ations have been attempted (see Bed ding & Kieldsen |2007| 
and lAerts et al.ll2008l for recent summaries), t Cet has the lowest 
mass. 

As a nearby bright star, t Cet has been intensively studied. 
A rotational period o f 34 days is suggested by sporadic peri- 
odicities in Ca II (Ba liunas et al ., 1996), but overall r Cet is a 
ver y inactive star w i th alrn ost no rotational modulation. This 
led lGrav & Baliunas (1 19941) to pr opose that t Cet is seen nearly 



pole-on, while lJudge et al.l(l2004 have suggested that it may be 
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in a phase analogous to the solar Maunder minimum. Its stabil- 
ity makes it a favoured targe t for testing the vel ocity stability of 
exoplanet programmes (e.g. But ler etakl [1996). Despite many 
velocity observations by different groups, no plan etary compan- 
ions have been reported (IWittenmver et al.L l2006l) . Direct imag- 
ing wit h the Hubble Space Tele scope also failed to detect acom- 
panion dSchroeder et aI.L l2000h . However. iGreaves et al. I (12004 
have imaged a debris disc around r Cet that has a dust mass at 
least an order of magnitude greater than in the Kuiper Belt. 

We note that t Cet is particularly suitable for an astero- 
seismic observing campaign because its radius has been de- 
termined interferometrically with an accuracy of 0.5%. The 
combination of interferometric and asteroseismic results has 
been applied to several o ther stars, as discussed in detail by 
Creevev et al. {20Q 3 and ICunha et all (|2007|) . As stressed by 
(Brown & Gillilan(jTl994 V for example, oscillation frequencies 
are most valuable for testing evolution theories when the other 
fundamental stellar properties are well-constrained, t Cet satis- 
fies this requirement as well as can be done for any single star. 



2. Velocity observations and guiding noise 

We were allocated six nights to observe r Cet on 2004 October 
2-7, using the HARPS spectrograph (High Accuracy Radial ve- 
locity Planet Searcher) on the 3.6-m telescope at the European 
Southern Observatory on La Silla in Chile. This spectrograph in- 
cludes a thorium emission lamp to provide a stable wavelength 
reference. 
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Time series of velocity measurements of t Cet. 
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Fig. 2. Velocity measurements of 6 Pav over 1 .5 hr at the start of night 3 
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We obtained 1962 spectra of r Cet, with a dead time of 31 s 
between exposures. For the first two nights we used an exposure 
time of 40 s (resulting in a Nyquist frequency of 7.04 mHz) but 
shortened this to 23 s (Nyquist frequency 9.26 mHz) for the re- 
mainder, in order to sample better the noise at high frequencies 
(see below). The velocities were processed using the method de- 
scribed by iBouchv et al.l (l200ll) and the resulting velocities are 
shown in Fig.[T] The fourth and fifth nights were mostly lost to 
poor weather 

For about 1.5 hr at the start of each night, when t Cet was 
inaccessible, we observed the star 6 Pav (HR7665; HD 190248; 
HIP99240; G6-8IV; V = 3.56). Small amounts of data on this 
star wer e also obtained with UVES at the VLT and UCLES at the 
AAT bv Kieldsen et al.k 2005). who found oscillations centred at 
2.3 mHz with peak amplitudes close to solar We obtained a total 
of 225 spectra of 6 Pav with HARPS, with exposure times of 
50 s (nights 1 and 2) and 23 s (nights 3 and 5). The velocities for 
night 3 (100 data points) are shown in Fig.|2] 

From the scatter in the velocities and the noise in the power 
spectra for both t Cet and 6 Pav, it was obvious that an unex- 
pected noise source was affecting the velocities. Figure [3] shows 
the power spectrum for night 3 for both stars. For r Cet (upper 
panel), there is a clear excess at 4 mHz, as expected for oscilla- 
tions in this star However, there is also a significant power ex- 
cess around 6 mHz. For 6 Pav (lower panel of Fig. the power 
centred at about 2.2 mHz is from oscillations (see lKieldsen et alj 
12005) but we again see additional power at 6 mHz. The effect 
of this instrumental noise is clearly visible by comparing the 
time ser ies of 6 P av in Fig. |2]of this paper with that in Fig. 2 
of .Kieldsen et all (2005). 

The spurious signal at 6 mHz was later traced to a periodic 
error in the telescope guiding system. Noise spikes at 3.1 and 
6.2 mHz have been reported in HARPS observations of other os- 
cillating s tars, namely 70 Oph (Carrier & Eggenberger, 2006), 
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Fig. 3. Power spectrum of velocity measurements for night 3 
only, for r Cet (9. 1 hr) and 6 Pav (1 .5 hr). 



In the case of t Cet, the 6 mHz noise is particularly problematic 
because it covers a fairly broad range of frequencies and because 
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Fig. 4. Power spectrum of t Cet for different subsets of the data. The left panels show the spectrum with the use of the raw weights, 
and the right panels show the spectrum after the weights had been adjusted to account for bad data points and for night-to-night 
variations in the noise level (see text). 



the stellar oscillations have very low amplitude. There does not 
seem to be a strong noise signal at 3 mHz in our data, although 
it is difficult to be certain. 



3. Data analysis and results 

Our analysis of the velocity data and of the extracted power 
spectrum for r Cet follows the method developed and used 
for a Cen A ( Butler et al., 2004; Bedding et al., 2004), a Cen B 
dKieldsen et al.> ,2005,) . v Ind (Bedding et al.> ,2006 ) and /? Hyi 



(iBedding et al.Ll2007l) . As usual, we have used the measurement 
uncertainties, cr,, as weights in calculating the power spectrum 
(according to w; - 1/crj). The main difference between rCet 
and other stars that we have analysed, apart from the problem of 
excess noise from the periodic guiding error, is the single-site na- 
ture of the observations. Because of this, we have not attempted 
to optimize the weights to reduce the sidelobes, in the way that 
we did for other stars. 
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3. 1 . Adjusting the weights 

As for previous stars, we adjusted the statistical weights to ac- 
count for bad data points and for night-to-night variations in the 
noise level. We did this by measuring the noise at frequencies 
where the oscillation signal and the long-term drifts are negli- 
gible. The first step was therefore to remove all power below 
800 yuHz (to avoid the slow drifts), as well as all power between 
3000 yuHz and 5500 /kHz (which is dominated by the oscilla- 
tions). This filtering was performed by the standard method of it- 
erative sine-wave fitting (sometimes known as 'pre-whitening'). 
In this method, the highest peak is identified in the region of the 
power spectrum that is to be removed, the corresponding sinu- 
soid is subtracted from the time series, the power spectrum is 
recomputed and the procedure is repeated until all the power is 
removed. 

Once this was done, the filtered time series for each night 
was examined for bad data points. These were identified as those 
deviating from the mean scatter by more than 4-cr, and were re- 
assigned lower statistical weights. We found that more than 10% 
of data points had to be significantly down-weighted. This frac- 
tion is much greater than for previous stars that we have observed 
and indicates the serious effects of the guiding errors. At the 
same time, we scaled the statistical weights on a night-to-night 
basis, in order to reflect the noise measured at high frequencies. 

The consequence of using the revised weights is a signifi- 
cant improvement in the signal-to-noise, as can be seen in Fig.|4] 
Comparison of the left and right panels in that figure shows that 
the adjustment of weights has removed essentially all the ex- 
cess noise at 6-7 mHz and also decreased the noise level in the 
range 1-3 mHz. The mean noise in those regions, measured in 
the amplitude spectrum of the whole data set, was reduced from 
6.0 to 4.0cm s '. Most of this reduction came from the down- 
weighting of bad data points, as described above. Also note that 
the strongest oscillation peaks are not as high in the combined 
data (top panel of Fig. |4]i as in the shorter subsets (middle and 
bottom panels). This reflects the finite lifetime of the modes (see 
SecO- 

3.2. The large separation 

The final power spectrum of r Cet is shown in the top-right panel 
of Fig.|4] There is a clear excess due to oscillations which is cen- 
tred at 4.5 mHz. The next step was to search for a regular series 
of peaks, as expected for p-mode oscillations, and to measure 
the large frequency separation, Av. We did this in two ways. The 
first was to smooth the power spectrum and then calculate the 
autocorrelation in the region of excess power, between 2.5 and 
6.0 mHz. This produced a clear peak at 169 /uHz. 

The second method, which we develop ed for the Kepler 
pipeline dChristensen-Dalsgaard et all l2007h . involved measur- 
ing the highest peak in the collapsed power spectrum for a range 
of values of the large separation. The collapsed power spectrum 
for a given value of Av is calculated by dividing the power spec- 
trum into intervals of length Av and summing these. The result 
is shown in Fig.|5] and again we see a peak at 169juHz. 

3.3. Osciilations frequencies and mode lifetime 

Mode frequencies for low-degree p-mode oscillations in main- 
sequence stars are well approximated by a regular series of 
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Fig. 5. Summed power of t Cet, using different values of the 
large separation. The dotted line marks the maximum, which oc- 
curs at a frequency separation of 169.2yuHz. 



peaks, with frequencies given by the following asymptotic re- 
lation: 

v„,, = Av(n +y + e)- 1(1 + 1)Dq. (1) 

Here n (the radial order) and I (the angular degree) are in- 
tegers, Av (the large separation) depends on the sound travel 
time across the whole star. Do is sensitive to the sound speed 
near the core and e is sensitive to the surface layers. See 
Christensen-Dalsgaard (2004) for a recent review of the theory 
of solar-like oscillations. 

We have extracted the individual oscillation frequencies for 
T Cet using the standard method of iterative sine-wave fitting 
down to S/N = 2.5. The single-site nature of the observations 
and the relatively low signal-to-noise ratio mean that this process 
is susceptible to one-cycle-per-day ambiguities (+11. 57 //Hz). 
On the other hand, we are helped greatly by the fact that both 
the large and small separations are much greater than the ex- 
pected mode linewidth, and so all modes are well separated. 
Furthermore, r Cet is an unevolved star and so we expect the 
oscillation frequencies to follow quite closely the asymptotic re- 
lation, without the presence of mixed modes. We have used this 
information to guide our selection of the correct peaks, but we 
stress that there is some uncertainty in the correct mode identifi- 
cation. 

The extracted frequencies are listed in Table[T] They are also 
shown in Fig. |6] in echelle format, where the frequencies are 
stacked in segments of length Av. As noted above, the mode 
identification is uncertain, and this is particularly true for the 
/ - and I = 2 modes above 5 mHz. 

A fit to these frequencies provides the various large and small 
separations, as listed in Table [2] For the definitions of these 
separations see iBedding & KieldsenI (l2003 b. for example. The 
separations generally vary with frequency and so the values in 
Table |2] are given at 4.3 mHz. The scatter of these frequencies 
about smooth ridges in the echelle diagram is about 1-2 /iHz, 
which indicates the uncertainties in the individual frequencies in 
the table. From this scatter we can estimate the mode lifetime, 
using the method described by Kjeldsen et al. (2005). We find 
a value of 1.7 + 0.5 d, which is sligh tly lower than the v alue of 
2.88 + 0.07 d measured for the Sun (iChaplin et al.Lll997h . 
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Table 1. Oscillation frequencies in t Cet (in fiHz) 



Table 2. Oscillation parameters for r Cet (see Sec. |3.3l l 
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Fig. 6. Echelle diagram of oscillation frequencies for r Cet, us- 
ing the frequencies listed in Table [U 



3.4. Stellar Parameters 

Detailed fitting of the oscillation frequencies of r Cet with the- 
oretical models is beyond the scope of this paper. However, we 
can use our results to determine the mean density of the star , 
via the empirical method described by ' Kjeldsen et al.l (l2008bl) . 
This method corrects the frequencies of stellar models for near- 
surface effects by making use of the fact that the offset between 
observations and models should tend to zero with decreasing fre- 
quency. The method involves fitting both Av and the absolute 
frequencies of the radial modes (i.e., those having degree I = 0). 



We applied the method to r Cet, using m odels computed with the 
Aarhu s stellar evolution code (ASTEC, IChristensen-DalsgaardI 
2008a') and the Aarhus adiabatic oscillation package (ADIPLS, 
Christensen-Dals2aard 2008b). The result is a mean density for 
T Cet of 2.21 + O.bl gcm-3 (0.45%). 

As mentioned in the Introduction, the angular diameter 
of T Cet ha s also been measured. This was first done by 
Pijpers et all fcOOS ^) with the VINCI instrument on the VLTI. 
They obtained an angular diameter, corrected for limb darken- 
ing, of 1.97 + 0.05 mas (2.5%), where the uncertainty was dom- 
inated by the uncertainty in the angular diameter of the cali- 
brator star. Subsequentlv. iDi Folco et alj 12004) used the VLTI 
with smaller calibrator stars to obtain an improved diameter of 
2.032 + 0.031 m as (1.5%). An e ven more accurate measurement 
was obtained bv lDi Folco et al.l ([2007) using the FLUOR instru- 
ment on the CHARA array, giving 2.015+0.01 1 mas (0.5%). The 
weighted mean of these measurements, which we adopt here, is 
2.022 + 0.010 mas (0.5%). Using the revised Hipparc os parallax 
for T Cet of 273.96 + 0.170mas (ivan Leeuwenll2067l) gives a ra- 
dius of 0.793 + 0.004 iJ0 (0.5%). Finally, combining this radius 
with our estimate from asteroseismology of the mean density 
gives a mass for r Cet of 0.783 ± 0.012 Mq (1.6%). 

For convenience, we also give an estimate of the luminosity 
of T Cet. The apparent visual magnitude ofV- 3.50 + 0.01, with 
the revised parallax, gives an absolute magnitude of My - 5.69 ± 
0.01. Using a bolo metric coiTection for rCet of -0.17 + 0.02 
jCasagrande et all |2006) and adopting an absolute bolometri c 
magnitude for the Sun of Mb^, © = 4.74 (iBessell et al.L Il998l) . 
we derive a luminosity for t Cet of 0.488 ± 0.010 L© (2.0%). 



3.5. Oscillation amplitudes 

We have determined the oscillation amplitude per mode from 
the smoothed power s pectrum, using the method described by 
iKieldsen et al.l (l2008al) . This produces a result that is indepen- 
dent of the stochastic nature of the excitation and damping. The 
result is shown in Fig. Q The peak of the envelope occurs at fre- 
quency v„iax - 4.49 mHz and the peak amplitude per mode (for 
radial oscillations) is 11.2 + 0.8cm/s, where the uncertainty is 
estimated using Eq. (3) of ' Kieldsen etaTI ( 12008 a'). For compar- 
ison, we also show in Fig.|7]the amplitude curves measured for 
the Sun (usin g stellar techniques) an d for a Cen B, both taken 
from Fig. 8 of IKieldsen et all (l2008al) . Note that thejuminosity 
and mass of a Cen B are, re s pectiv ely, 0.51 Lq (lYildizl |2007|) 



ss pectiy i 

and 0.93 Mq (Pourbai x et all. l2002h . With this in mind, we see 
that the relative amplitudes of the three stars in Fig.Qare in rea- 
sonable agreement with th e L/M s cahng relation prop osed by 
IKieldsen & Bedding! (I1995L see also lSamadi eTani2007l) 
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Fig. 7. Smoothed amplitude curve for t Cet. For compar- 
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4. Conclusions 

We have used HARPS to measure oscillations in the low-mass 
star T Cet. Although the data were compromised by instrumen- 
tal noise, we have been able to extract the main features of the 
oscillations. We found t Cet to oscillate with an amplitude that 
is about half that of the Sun, and with a mode lifetime that 
is slightly smaller than solar The large frequency separation 
is 169 //Hz, and we have identified modes with degrees 0, 1, 
2 and 3. It is important to stress that, given the relatively low 
signal-to-noise ratio and the single-site nature of the observa- 
tions, there is some uncertainty in the correct mode identifica- 
tion. 

We used the frequencies of the radial modes to estimate the 
mean density of the star to an accuracy of 0.45%, from which 
we derived a mass of 0.783 + 0.012 M© (1.6%). More detailed 
modelling of the oscillation frequencies will be the subject of a 
future paper. 
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